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Abstract 
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1 Introduction 



In the limit of vanishing hght quark masses m„, and the QCD Lagrangian exhibits a 
symmetry under chiral SU{3)l x SU{3)ji transformations which is broken down spontaneously 
to SU{3)v giving rise to the eight pseudoscalar Goldstone bosons {7!-,K,ri). The interactions 
at low energies amongst the Goldstone bosons can be treated in a perturbative fashion within 
the framework of chiral perturbation theory (ChPT). The axial U{1) anomaly of the QCD 
Lagrangian, on the other hand, prevents the corresponding pseudoscalar singlet from being a 
Goldstone boson. Indeed, the lightest candidate would be the r]' with a mass of 958 MeV which 
is considerably heavier than the octet states. Therefore, in conventional ChPT the t]' is not 
included explicitly, although its effects are hidden in a contribution to a coupling coefficient of 
the higher order Lagrangian, a so-called low-energy constant (LEG). 

In order to describe the physical processes involving the t]', such as the hadronic decay modes 
7]' — i> T]nn, TTTTTT, the anomalous decays rj' — >■ 77, 7r+7r~7, as well as photoproduction of the t]' 
on the nucleons and r]' production in pp collisions, it is necessary to include the r]' explicitly 
within the effective theory. The inclusion of the t]' can be achieved either by introducing the 
topological charge density as a background field in the effective Lagrangian, see e.g. [|l], |[, or by 
extending the chiral SU{3)l x SU{3)n symmetry to U{3)l x U{3)r |^, ^, 0. (The equivalence 
of both schemes has been proved in 0.) 

Due to the large t]' mass, however, the convergence of the chiral expansion is questionable. 
As a matter of fact, the appearance of a mass that does not vanish in the chiral limit, such as 
the mass of the rj' will spoil the chiral counting scheme so that higher loops will still contribute 
to lower orders. This can be prevented by imposing large Nc counting rules within the effective 
theory in addition to the chiral counting scheme. In the large N^. limit the axial anomaly 
vanishes and the t]' converts into a Goldstone boson. The properties of the theory may then 
be analyzed in a triple expansion in powers of small momenta, quark masses and l/N^, see 
e.g. IQ, ^ 1^ ^. In particular, m^/ is treated as a small quantity. Phenomenologically, this 
is not the case and we will therefore treat the rj' as a massive state leading to a situation 
analogous to GhPT with baryons. Recently, a new regularization scheme - the so-called infrared 
regularization - has been proposed which maintains Lorentz and chiral invariance explicitly at 
all stages of the calculation while providing a systematic counting scheme for the evaluation of 
the chiral loops 0. Infrared regularization can be applied in the presence of any massive state 
and has been employed in U{3) GhPT in [0, [Tl[| . 

In the present work, we apply this approach to the dominant decay of the 77': 77' — t^tttt, 
including one-loop corrections. This decay has been the subject of many investigations using 
chiral effective Lagrangians, see e.g. |T3|, |T^, |T5|, however one-loop corrections have not 

been included so far. Of particular interest is the recent work ||15[ where a systematic study up 
to next-to-leading order within the framework of large GhPT has been undertaken. One-loop 
corrections have not been included in since they start at third order in the combined chiral 
and 1/Nc expansion, and although the results at tree level were promising, no statement about 
the convergence of the expansion could be made. If, on the other hand, large counting rules 
are not employed, loop corrections start already at next-to-leading order. Our investigation 
will shed some light on the convergence issue of infrared regularized U{3) GhPT. It includes for 
the first time unitarity corrections whereas in |jlO|, |ll| tadpoles were the only loop corrections. 
The present calculation will be a crucial test for the convergence of the chiral series within this 
approach. 
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It has been claimed that this decay can be described in a tree-level model via the exchange 
of the scalar mesons (t(560), /o(980) and ao(980) which are combined together with k(900) 
into a nonet [|T^ (see also [0, |T^, ^ |^)- The authors find that the exchange of the 
scalar resonance ao(980) dominates. It will be interesting to see if we are able to reproduce 
the experimental data reasonably well without the explicit inclusion of resonances such as the 
ao(980) although their effects might be hidden in some LECs. 

In the next section we present the formalism for the rj rjini decay as well as the pertinent 
Lagrangian up to fourth chiral order and contributions from tree diagrams. The appearing chiral 
loop corrections are discussed in Sec. |^ and Sec. § contains the comparison with experiment. 
We summarize our findings and conclusions in Sec. ^. 



2 Formalism and leading contributions 



The process rj' rjiT-K is the dominant decay mode of the rj'. The measured rates are ||22 

r(V ^ r/TT+TT") = (90 ± 8) keV 



with the ratio given by 



T{r]' r77r°7r°) = (42 ± 4) keV (1) 
r(rj' rj7r~^7i~ 



2.1 ±0.4 (2) 



which is consistent with the ratio of 2 in the isospin limit. Isospin violating corrections to the 
decay rates seem to be small and will therefore be neglected throughout this work. Taking this 
into account the amplitudes of the two decay channels are equal 

^ri'-^riTT+TT- (s, t, u) = A^i _^^T^o^O (s, t, u) (3) 

where we introduced the Mandelstam variables 

S = iPn' - Pr,)^ t = {prj' - U = (p^/ - p^a)^ (4) 

with being the four-momentum of the meson if. In order to calculate the amplitude, we 
restrict ourselves to the effective Lagrangian up to next-to-leading order, i.e. to fourth order in 
the chiral expansion. The construction of the pertinent Lagrangian has been given elsewhere, 
see e.g. [^, |^, ||, |T0|, and we sketch the results briefly. Note that l/N^ counting rules are 
not necessary to set up the Lagrangian including the t]'. The effective Lagrangian for the 
pseudoscalar meson nonet {n, K,ris,riQ) reads up to second order in the derivative expansion 

C^'^'^ = -Vo{Vo) + EMVo)Ol'^ (5) 

with the relevant second order operators 

Of) = -(C^C^), Of = (M), 0f = (ziV), 0? = {Cn{C,). (6) 

We used the abbreviations = Wd^U, M = Wx + X^U, N = ^ X^U, and f/ is a unitary 
3x3 matrix containing the meson nonet. Its dependence on rjg and rjo is given by 

U = exp(diag(l, 1, -2) ■ tr]s/V3f + zv^r^o/Vs/ + •••)• (7) 
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The expression (...) denotes the trace in flavor space, / is the pseudoscalar decay constant in 
the chiral hmit and the quark mass matrix Ai = dia.g{'mu,rnd,ms) enters in the combination 
X = 2BM. with B = — (0|gg|0)//^ being the order parameter of the spontaneous symmetry 
violation. In the present investigation we work in the isospin limit m„ = 771^ = m and, 
therefore, only 77-77' mixing occurs. We have neglected terms which do not contribute to the 
process considered here, such as those containing external sources. 

The coefficients Vi are functions of the singlet field and are expanded in terms of 770 

„.M=„<n„<')^+„f)^+„f)5|+... (8) 

At a given order of derivatives of the meson fields U and insertions of the quark mass matrix Ai 
one obtains an infinite string of increasing powers of rjo with couplings v^''^ which are not fixed by 
chiral symmetry and have to be determined by phenomenologically. Parity conservation implies 
that the Vi are all even functions of 7/0 except fs, which is odd. The correct normalization for 
the quadratic terms of the mesons is obtained by vi{0) = V2{0) = vi{0) - 3v^{0) = \p. 
The next-to-leading order Lagrangian reads 

C^'^ = E^/3^iV0)O, (9) 

with contributing fourth order operators 

O = {C^^C^C.O, O 1 = (C'^C^)(C"^a), 

On = -{C>'){C^C-C,), Oi4 = -{C^^){C,){C^C,), 

Oi5 = -(c'')(c"')(c^a), = (c^)(c^)(co(a), 

O 4 = -(C^C^)(M), O 5 = -(W^M), 

Oi7 = {C^){C,){M), O18 = -{C^){C,M), (10) 

O21 = {C^^C^iN), O22 = {C^^C^){tN), 

O23 = {C^^){C,iN), O24 = {C^^){C,){iN), 

0, = {M){M), 7 = {N){N), 

O s = ^{MM + NN), Oi2 = l{MM-NN), 

02, = {tMN), 026 = (M)(ziV), 

where we kept the notation from ||^ . The coefficients jSi are expanded in tjq in the same manner 
as the Vi in (||); they are even (odd) functions of r/o for an even (odd) number of operators 
N. It turns out to be more convenient to include the /3o term, although there is a Cayley- 
Hamilton matrix identity that enables one to remove this term leading to modified coefficients 
Pi, i = 1,2, 3, 13, 14, 15, 16 0. Here we do not make use of the Cayley-Hamilton identity and 
keep all couplings in order to present the most general expressions in terms of these parameters. 
One can then drop one of the Pi involved in the Cayley-Hamilton identity at any stage of the 
calculation. 

Having presented the effective Lagrangian up to fourth chiral order we now proceed by 
calculating tree diagrams contributing to this decay. The contribution from C^'^^ reads 

_ AV2vi'^ml 

A. - (11) 
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with = 138 MeV being the average pion mass, /^r = 92.4 MeV the pion decay constant and 
the abbreviation 



V? = If. - IVQv? . (12) 

Inserting the value ^2^'* = 2.67 x 10^^ GeV^, obtained from a fit to the pseudoscalar decay 
constants [|l^, leads to a decay width of F = 1.9 keV which is about 20 times smaller than the 



experimental value. This has already been pointed out in [O, O, |T^, where an even smaller 



result is obtained due to the omission of ^3^'*. Note that we have replaced the pseudoscalar 
decay constant / appearing in the Lagrangian by the physical parameter f-j^ which coincides 
with / in the chiral limit. The difference in the decay amplitude is of higher order and is 
included in the higher chiral order contributions, see below. We proceed in a similar way by 
substituting the physical pion and kaon masses for the corresponding masses at leading order 
and writing the differences into the higher order contributions. An important feature of ([TT|) 
is that 77-77' mixing does not show up at this order and, as a matter of fact, the amplitude 
calculated in ( pT]) corresponds to the decay 770 — t^stttt. Within the used scheme rj-rj' mixing is 
of second chiral order |jll| and the corrections due to mixing on the decay amplitude ([Tl|) will 



be of fourth chiral order. 

Apparently, tree level diagrams from the Lagrangian 

£(2) 

are not sufficient to obtain reason- 
able agreement with experiment. One must therefore go beyond the leading order Lagrangian 
and take into account tree diagrams from the Lagrangian C^^^ and loops with vertices from 
The loop contributions will be discussed in detail in the next section and we will start here by 
evaluating tree diagrams from C^^\ The complete tree level amplitude A4 from C^^^ is given in 
App. |B[ Comparing the tree level contributions one makes the following crucial observation: 
although, strictly speaking, of fourth chiral order in the expansion in the derivatives and meson 
masses, the terms with the LEG combinations /3o,3,i3, /55,i8, /^l^^ and ^^^^ (see App. for the 
definitions) contribute already at lower orders if one makes use of p"^^, = m^, ~ 0{1). More 
generally, we will consider the Mandelstam variables s, t and u as zeroth chiral order. This 
counting is motivated by the fact that in the chiral limit s,t and u range from to m^,. From 

that point of view, it is more convenient to summarize the terms ^2^^ /3o,3,i3; /^5,i8; /^l^"* and j3^^ 
in 

4a/2^(i) 2 2 2 

+ ^/9o,3,i3(3^ + e + u'- 2ml - K - O 

16v^(7n|-m2)/3f) , 2 N 2 . N 

4^(2) 

J TT 

and to distinguish them from the remaining terms in A4. The other contributions from the 
fourth order Lagrangian which have not been included in (O) proportional to m'^ 

with if = n, K,r] and therefore contribute to the decay amplitude at fourth chiral order. This 
clearly separates the terms given in (0) from the remaining tree level contributions in App. 
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The last two terms are due to rj-rj' mixing: Asstttt /Jq ft denotes the contributions of the 

>(2) 



amphtude Asstttt in App. |B| proportional to /9q°^ + and and -RS, R^q^ describe the 



mixing of the rj and rj' fields at one-loop order ||rT| with 

< = |v^«--:)4V(/>o), 

?(2) _ / 10_2^(0) _ .2 MO) 



?{2) 
•-Or)' 



<' = -4m^/54,5,i7,i8//', (14) 



where m^ = 2v^^ /f^ and mg = + The chiral logarithm Ak is given in (0). We 

observe that the contribution proportional to /3o,3,i3 is the only term in that does not vanish 
for zero quark masses; for reasons to be explained below it is expected to dominate. 

In order to obtain a rough estimate for the tree diagram contributions, we will borrow 
the values for the coefficients (3^^^ from conventional SU{3) ChPT pSf. There is, however, 
one complication: in SU{3) ChPT the operator Oq is usually eliminated by making use of a 
Cayley-Hamilon identity which changes the l/N^ hierarchy of some of the LECs Lj We 
prefer to keep it explicitly but must transform the values of the Lj accordingly. An estimate 
of the coupling of Oq that respects the 1/Nc hierarchy is given by QCD bosonization models 



24, m 



^° 1927r2 ' ^ ^ 

In order to introduce the /Jg term within the SU (3) framework one must use the transformed 
couplings p^"'^^^ as follows 

/3f = Li - i/?o'''''\ Pr'' = L2- Po""''', Pr'' = L^ + Pr'' = U for I > 4. 

(16) 

Note that P^^^^^ and /^l*^^^^ are both considerably closer to zero than Li and L2 in agreement 
with large A*",, rules. The values pf''^ can be taken directly from the transformed values of SU{3) 
ChPT with the exception of l3^^\ The value of jSj^*^^^ is entirely saturated by the exchange of 



an rj' p8[. Starting with = in the extended framework and integrating out the rj' field 
produces the value of P^^^^^ |Tl|. Hence we set 



Having fixed the LECs which already appear in conventional SU (3) ChPT we can now proceed 
in estimating the remaining pf^^ with i > 13. To this end, we notice that these new unknown 
couplings enter only in combinations with the known pI^\ i = 0, . . . ,8, e.g. P^ and P^ appear 
in /5o,3,i3 and P^^is, respectively, where they represent OZI violating corrections (for P13 this is 
only true if the Cayley-Hamilton identity is used as described above). Therefore, we make the 
rough estimate by neglecting the unknown OZI suppressed couplings within these parameter 
combinations and keeping only the known parameters. For the two cases mentioned this yields 
/?o,3,i3 ~ f^o'^ + /^s^'* and P^^ig ~ P^^ . Furthermore, P^^ and pf^ are OZI violating corrections 
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to (3^^ and we use the estimate Pa^^^u^is ~ <i(3f'^ + (3^^ This rough estimate for the parameter 
combinations is motivated by conventional SU (3) ChPT, where a decent estimate of the LECs 
can be obtained via resonance exchange if resonance couphngs are used which obey the OZI 



rule [|28|. Contributions from OZI suppressed couplings to the operators in SU{3) ChPT seem 
to be negligible and we will assume the same also in U{3) ChPT. 

Finally, the parameters /32i, /?22 and P25, P26 from parity-odd operators enter the calculation. 
The coupling P21 appears together with P^, while P22 appears with /34. Both P21 and P22 are 
suppressed by one power of 1/iVc in comparison with jS^ and P^, respectively. The same applies 
for P25 {P26) which combines with /?§ (/Jg + Pj)- We will also drop these 1/Nc corrections from 
parity-odd operators, hence generalizing our approximation for the OZI violating contributions. 
In particular, and P^^^ are approximated by P4 and P5, respectively. Making use of the 
1/Nc hierarchy for the LECs, we are only left with known parameters pf'^ of the fourth order 
Lagrangian. 

The same estimate of dropping 1/Nc corrections in parameter combinations can also be 
applied to coefficients from the lowest order Lagrangian. Keeping = 2.67 x lO'^ GeV^ 



(2) 

fixed at its value determined from a fit to the pseudoscalar decay constants ^2 enters 
only in the combination where it is 1/Nc suppressed. It will thus be neglected, yielding 
the estimate v^^ ~ v^^ — ~ 3.2 x 10^'^ GeV^. We also omit the coefficient v^^ which 

should have a small value since it is 1/Nc suppressed with respect to v^^ = P /A. Indeed, small 
variations of v\ within a realistic range have shown that this contribution has almost no effect 
on the results. Dropping 1/Nc suppressed corrections simplifies the parameter combinations 
considerably and reduces them to the contribution of the leading Nc piece. Afterwards we 
are mainly left with known parameters from conventional SU{?)) ChPT. Of course, one is free 
to keep the 1/Nc corrections, but our results will lose most of their predictive power and it 
will be relatively easy to accommodate the experimental data by finetuning these LECs. The 
neglection of such couplings provides a more stringent test of our approach. Estimating the 
LEC combinations in the amplitude by the contributions from the leading Nc pieces should 
not be confused with imposing large Nc counting rules in ChPT. We do not make use of large 
Nc rules in the calculation of the amplitude which is obtained by applying infrared regularized 
C/(3) ChPT. 



Values for the LECs in S'f/(3) ChPT were given in p3[. They correspond to the values 
= 1.6 ± 0.0, pf^ = -0.4 ± 0.3, ^ = -0.2 ± 0.3, 4°^ = -0.3 ± LI, /?f ^ = -0.3 ± 0.5, 
/3f ^ = 1.4±0.5, = -0.2±0.3, = 0.0±0.2 and p^ = 0.9±0.3 (in units of lO^^) in U{3) 
ChPT when applying the procedure described above. All the other LECs are approximately 
zero and will be neglected. The three combinations of LECs in App. have the values 
/?o,3,i3 = 1.2 ± LI, /35,i8 = L4 ± 0.5 and /54,5,i7,i8 = 0.5 ± L6. 

We will give amplitudes at a central point in phase space, the Mandelstam variables being 
given by their phase space expectation values 

So = 0.12GeV^ to =^^0 = 0.567 GeVl (18) 

The phase space of the decay is rather small, so that fiuctuations of the amplitude have only 
little impact on the decay width and incidentally the squared amplitude at this point is a good 
approximation for the decay width in units of keV. 

In order to obtain numerical results, we first note that the amplitude is dominated by the 
/9o,3,i3 term. The huge uncertainty in the value of /5o,3,i3, which is as large as the value itself. 
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makes a numerical prediction of the decay width impossible, it only allows for an upper bound 
of roughly 200 keV. Hence we will use the experimental decay width and the slope parameter 
a (see below) to fit some of the LECs and reduce the uncertainty of /9o,3,i3. We find that 
the values + = 1.1, /5f^ = -0.7 and = 1.9 (in units of lO'^), which lie inside 
their uncertainties given above, reproduce the data well. Alternatively, one could employ the 
central values but finetune the OZI violating corrections which we have neglected initially. The 
contributions from l/N^ suppressed parameters are then within the given error ranges for the 
known LECs and our initial assumption of omitting the 1/Nc corrections seems to be justified. 
Our estimate for the amplitude reads then 

^Lo(so,to,^o) = -6.4. (19) 

For comparison the decay width integrated over phase space is 

|r(r/' ^ r]TT+TT-) = r(V ^ r/7r°7r°) = 42.2 keV. (20) 

which coincides numerically with \ Ai^o{so,to, uo)P ~ (6-4)^ ~ 40.2. Breaking down into its 
constituents one obtains 

Ao(so, to, uo) = 1.3 - 1.7 - 7.8 + 1.3 + 0.5 + 0.0 (21) 

with the 52^\/35,i8,/5o,3,i3,^f\/35^^ and aJJ^stttt terms in order. Apparently, the /3o,3,i3 term 
dominates the amplitude being an order of magnitude larger than v!^^ and /Ss^is- In previous 
work it was realized that the only contribution to the decay from the lowest order Lagrangian 
£(2) jg given by v^^ yielding a width much smaller than the measured one []T2|, [T^. This 



dilemma was usually cured by adding momentum dependent terms of the type Po . . . 1 12, 14 



within the framework of large Nc ChPT the /3o term turned out to be dominating although it 



is a next-to- leading order effect [15 



After fitting the value of /3o,3,i3 to the decay width we may reduce its uncertainty to 

/3o,3,i3 = (l.l±0.4)xlO-l (22) 

The error accounts for the uncertainties of both the decay widths and the other LECs. 

The dominance of momentum dependent terms can be easily understood by the requirement 
of Adler zeros. Assuming analyticity in the Mandelstam variables s, t and u (which is certainly 
the case for tree diagrams), the amphtude can be expanded as 

A{t, u) = a + b{t + u) + c{t + uf + d{t - uf + higher orders in t, u (23) 

where we made use of the symmetry under the exchange of the two pions. Since for 0, i.e. 
s = m'^,,t = ml,u = rri^, the amplitude vanishes, a must be proportional to m^. For physical 
values of the Mandelstam variables, on the other hand, one has the constraint t+u > m^,. From 
the ratio (m^/m^')^ ~ 0.02 it becomes obvious that in the physically allowed region momentum 
dependent terms might be important and even dominate the leading order contribution a. 

One could now argue that higher order terms with an increasing number of derivatives on 
the singlet field r/o will lead to contributions of the type (t + u)"' and that consequently the series 
in (p3D will not converge. To this end, we would like to take a closer look at terms with more 
derivatives on the singlet field rjQ. Neglecting small corrections from rj-r]' mixing, possible tree 
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level contributions have the generic form d''rjod'"rjsd'^7cd"'7T where we omitted Lorentz indices for 
brevity. Derivatives acting on one field cannot have the same Lorentz index since such terms 
are eliminated by using the equation of motion for the mesons. This means that derivatives on 
r]Q must be contracted with derivatives on a Goldstone boson which are suppressed by one order 
in the Goldstone boson masses. Hence, we expect the series to converge, albeit slowly. In order 
to clarify this issue the inclusion of terms from the sixth order chiral Lagrangian, in particular 
those with six derivatives, would be very helpful. However, at that order two-loop diagrams 
contribute together with a proliferation of many new counterterms. Such an analysis is far 
beyond the scope of the present work. We will therefore restrict ourselves to the Lagrangian up 
to fourth chiral order, assuming that corrections from higher orders in the derivative expansion 
are well-behaved. 



3 Chiral loops 

To our knowledge, the calculation of chiral loops has been omitted in previous work. This is 
due to the fact that the rj' propagation within a loop introduces a hadronic mass scale, m^/, and 
spoils the chiral counting scheme if dimensional regularization is employed for the evaluation 
of the integrals. This can be prevented using large ChPT and a first step towards this 
direction has been undertaken by calculating the next-to-leading order tree diagram corrections 



ISf . Within this scenario loop diagrams start contributing at third order in the combined chiral 
and 1/Nc expansion. Although the tree level results were quite promising, it remains to be seen 
whether loop corrections destroy the agreement with experimental data, especially, since the t]' 
mass is treated as a small quantity which is, in our opinion, phenomenologically not justified. 
As mentioned before, we treat the t]' as a massive state and thus the situation is similar to 
the case with baryons in conventional ChPT. In order to establish a chiral counting scheme 
even in the presence of massive states, a new regularization scheme has been proposed in |^ 



following the ideas of [EO, |3D|. The so-called infrared regularization keeps Lorentz and chiral 



invariance explicit at all stages of the calculation and is thus perfectly suited for U{3) ChPT. 
The underlying idea of this regularization scheme is to extract the infrared singularities of the 
loop graphs in a relativistically invariant fashion. The chiral expansion of the remainder of 
the integral is an ordinary Taylor series in the momenta and Goldstone boson masses and is 
absorbed into the parameters of the counterterms which are renormalized accordingly. This 
means that the finite pieces of the infrared part of the integral are kept whereas the remainder 
is omitted. 

Infrared regularization has already been applied within U{3) ChPT in |ll| where the 
resulting chiral series appeared to converge. In both investigations only tadpoles contributed 
at one-loop order. In order to make a more profound statement about the convergence of the 
chiral expansion , it is necessary to consider unitarity corrections. Indeed, for the decay process 
considered in this paper both tadpoles and one-loop diagrams of the self-energy type contribute, 
see Fig. ^ and we will calculate the leading chiral loops, i.e. the leading chiral corrections from 
loops with the vertices (^5^18 and /9o,3,i3- 

The results for the fundamental integrals have already been presented in which we 
briefly discuss here restricting ourselves to the calculation of the finite pieces and neglecting 
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Figure 1: s, t, u and tadpole loop diagrams contributing to the decay. Solid lines denote the 
pseudoscalar mesons. 



the divergences. The finite infrared parts of the tadpoles are readily evaluated 

dH i 



2 2 



dH i 



I {27rY P - ml, + ie 



= 



(24) 



where the subscript / denotes the infrared portion of the integral and ^ is the scale introduced 
in infrared regularization. The vanishing of the rj' tadpole is due to the lack of a low-energy 
scale within the integral, such as the Goldstone boson masses. 
The fundamental scalar integral with two propagators reads 



dH 



J {27rY {{p - If -ml + ie) {P - ml + ie) 



(25) 



If both mesons inside the loop are Goldstone bosons the infrared part is obtained by regularizing 
the integral dimensionally. Expressed in terms of elementary functions the regularized integral 
equals 



167r2 



- 1 + ln 



mamb 



+ 



2^/A 



artanh 



A 



(ma + m^y — p^ 
A = {{ma - m^f - p^) {{ma + m^f - p^) . 



(26) 



Outside their real- valued domains (x > 0) the involved functions may be replaced according to 
\/—x — i^/x^ artanh = i arctanx, artanh(— 1 — x) = |i7r -|- arcoth(— 1 — x). (27) 
For a loop with an 77' and a Goldstone boson the pertinent result reads {a — tt, K,r]) 

1 a{n + a) 



Iv'a{P ) 



167r2 1 + 2aVt + 



(Ina'-l) 



+ 



1 , n + a + ^/W^ 
m 



167r2 1 + 2q;Q + q;2 + q; - 



(28) 
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where a = ma/rrir^i and Vt = {jp' — rri^, — m^)/ (2m^/ma) is counted as a zeroth order quantity, 
Vt = 0(5°), and we have omitted the divergent pieces. The result ( ]28|) is valid for > 1; for 
< 1 analytic continuation yields 

h'aip ) = n ^ — 7z — ^(Ina^ - 1) 

+ — ^ ^ arccos = (29) 

which is the result presented in 0. 

The chiral expansion of this expression has been shown to converge in the interval 

< ^ < • (30) 



2a 2a 

Since we do not consider isospin violation, a loop with an r\ and rj can only occur in the s 
channel. In this case, f2 is constrained by 

1 + a2 

-l>f]>-1.09> — ^-1.16. (31) 

2a 

In the t and u channels the r( propagates only in combination with a pion and again the 
constraint for the limits is fulfilled 

1 + a2 

-1 > ^] > -1.77 > — ^ -3.54. (32) 

- 2a ^ ' 

Hence, in both possible cases with exactly one r( inside the loop the chiral series converges and 
we can expand the expression (^) in a at fixed fi, obtaining /r,'a(p^) = 0{a). In comparison 
with the pure Goldstone boson loop which is of chiral order 0{\\ the case with one r/' propagator 
starts contributing at one chiral order higher. In other words, the leading chiral contributions 
from the loops are the ones with two Goldstone bosons, while those with an Vj and a Goldstone 
boson are suppressed by one chiral order. It is thus beyond our working precision of taking 
only the leading chiral loops into account and we can safely neglect these contributions. 

Finally, in the s channel two v( mesons can propagate inside the loop. This integral was 
shown in |^ to contain no infrared singularities and we will omit these loops as well. In order 
to calculate the leading loop corrections, it is sufficient to consider the pure Goldstone boson 
loops, while dropping those with the propagation of an Vj . To the order we are working the 
77' field can therefore be regarded as an external field which does not appear in loops. This 
argument can be turned around by treating the Vj field as an external field right from the 
beginning and considering Goldstone boson loops only. In the present calculation effects of the 
rl within the loops appear at sub-leading order whereas they are omitted completely in the 
latter approach. 

In App. ^ we have presented the amplitudes including the leading non-analytic contribu- 
tions in infrared regularization. The terms and are the tree diagrams from the corre- 
sponding Lagrangians ^C-^) and ^C-^^ Az stems from the wave function renormalization of the 
octet fields, A^/ is due to the higher order corrections of the meson masses and the pion de- 
cay constant and Ai^^Ag^t.u are the tadpole and the unitarity corrections in the s, t, w-channel, 
respectively. We have presented the amplitudes for the bare unmixed fields r^o and 773, e.g. 
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A-ssiT-K = ^sstttt + ^sstttt is the expression for two external ris fields. In order to obtain the 
amplitude for the physical fields rj and rj', one must include rj-r]' mixing and wave-function 
renormalization of the tiq. It would then be sufficient to calculate Aqstttt up to fourth chiral 

(2) (2) 

order and both Asstttt and Aootttt only at lowest order A^^'^^, Aqq^^ and take rj-rj' mixing into 
account. However, this result turns out to be scale dependent both at second and fourth chiral 
order due to s,t,u ~ 0{1). The scale dependence at second order is removed by calculating 
the complete one-loop amplitude A^gg^^ as given in App. |B| and applying rj-r]' mixing at second 
order. The residual scale dependence is then of fourth chiral order and proportional to the 
mixing angle. It could in principle be cancelled by including, among other contributions, r]-r]' 
mixing at fourth order which is a two-loop effect and beyond the accuracy of this calculation. 
We will therefore content ourselves by removing the scale dependence at second chiral order 
and leave the scale dependence at fourth order which is numerically less significant. Note that 
the amplitude Aootttt does not exhibit any scale dependence up to fourth chiral order and it 

(2) 

suffices to include only the lowest order contribution AgQ^^g- The complete amplitude we are 
using is given by 

^NLO = (l + -Ror?^' ) ^OStttt + Rgiji {AsSnn — Rg^ AsSmr,2) + -Ror?^^007r7r,2- (33) 

The factor uf^, accounts for the Z-factor of the rj' field, Uf^, , R^^^ are due to r]-ri' mixing, and 

(2) 

i?g^ corrects for the Z-factor of the ris, since it was applied to both r/g fields in Asstttt. Using 
the rough estimates for the parameters (3i as explained above we obtain at the central point 
in phase space, {so,tQ,uo), for /i = m^/ AnlqI/^ = "^j?') = —5.1 + 0.22z and for /i = nip the 
result reads A^i^oif^ — ^p) — ~4.4 + 0.22i. The values of the non-zero LECs /?f ^ were shifted 
according to their renormalization constants to accomodate for a modified scale /i. 

However, a fit to the experimental value of the decay width can easily be accommodated 
by varying the LECs within their phenomenological ranges, e.g. /3o,3,i3 can be modified while 
keeping the other couplings fixed. The values of (/9o,3,i3 = 1-4 at = nip) or (/9o,3,i3 = 1.3 at /i = 
rriri') lead both to Anlq = —6.6+0.222. Such a fit is not unique since it could also be accomodated 
by varying some of the other LECs within their penomeno logical ranges and/or finetuning the 
OZI suppressed corrections and we will refrain from performing one. Furthermore, some of the 
higher order corrections (see below) tend to increase the magnitude of the amplitude cancelling 
partially the effect of A^^^q. 

Clearly, the loop corrections using the couplings from the leading Lagrangian 

£(2) 

are under 

control with a scale dependence which is considerably smaller than the uncertainty of the 
parameters see (p^Qj). 

3.1 Loops from the fourth order lagrangian 

In the preceding section, we have included loop contributions with vertices from the lowest order 
Lagrangian C^'^\ However, certain terms from the Lagrangian 

£(4) 

appear to spoil the strict 

chiral counting scheme due to s, t, m ~ C'(l) and contribute at lower chiral orders as one would 
naively expect. In fact, under reasonable assumptions for the LECs the /3o,3,i3 term dominates 
the other tree level contributions by an order of magnitude. It is therefore a legitimate question 
to ask, whether loop diagrams with vertices from in particular /9o,3,i3) are numerically under 
control and do not upset the convergence of the amplitude expanded in terms of derivatives 
and Goldstone boson masses. In this section we will clarify the issue by calculating one-loop 
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diag. 




ap 






V 


+1.4 + 0.02 


+0.00 + O.OO'i 


+1.4 + O.Oi 


+0.00 + 0.00i 


Po 


-8.6 + O.Oz 


-0.19 + O.OOi 


-8.6 + O.Oi 


-0.19 + O.OOi 


Pa 


+ 1.4 + 0.0i 


+0.12 + 0.00i 


+1.4 + O.Oi 


+0.12 + 0.00i 




-1.1 + 0.0z 


-0.00 + O.OOi 


-1.1 + O.Oi 


+0.00 + 0.00i 


v—v 


+0.9 + 0.2i 


-0.03 + 0.03Z 


+1.7 + 0.2i 


-0.07 + 0.03i 


I3q-v 


-0.1 - l.Oi 


-0.09 - 0.14z 


+0.0- l.Oi 


-0.11 - 0.14i 


(3a-v 


-0.1 + 0.2i 


-0.09 - 0.20i 


-0.6 + 0.2i 


+0.25 -0.20i 


P5-V 


-0.2 - 0.2i 


+0.02 - 0.02i 


+0.6 - 0.2i 


+0.06 -0.02i 


Po-Po 


-0.0 - 0.2i 


+0.01 + 0.00i 


+0.1 - 0.2i 


+0.02 - O.OOi 


Po-Pa 


-0.0 + O.li 


-0.51 - 0.30i 


+0.1 + O.li 


-0.76 - 0.30i 




+0.2 - O.Oi 


-0.05-O.Oli 


-0.2 - O.Oi 


-0.08-O.Oli 


Pa-Pa 


-0.1 -O.Oi 


-0.08 + O.OOi 


-0.0 - O.Oi 


-0.19 + 0.09i 


Pa-P, 


+ 1.1 + O.Oi 


+0.18-0.06i 


+1.8 + O.Oi 


+0.60 - 0.05i 


P5-P5 


-1.4 - O.Oi 


+0.03 - O.OOi 


-2.9 - O.Oi 


+0.06 + O.OOi 



Table 1: Amplitudes from higher loops, "a-6" denotes a loop between an a and a h vertex, "w" 
denotes any vertex from the second order Lagrangian. For comparison we give the tree level 
amplitudes "a" of the considered vertices a. A gives the amplitude at the central point in phase 
space so,to5^o ( |l8l) and a = {— 1/^/42) {OA/ dy) corresponds roughly to the slope parameter a 
of (|36|). The loop integrals are evaluated at the scales /x = rrip and /i = m^/. 

diagrams with /3q°\ P^^ and pf^ vertices, which were part of the /^s^is, /3o,3,i35 Pa^^ and ^5^^ terms 
we have included in the leading order amplitude Ai^o, (P^- Strictly speaking, they are of fourth 
chiral order in the derivative/meson mass expansion of the effective Lagrangian. This implies 
that the one-loop diagrams with these vertices have eight powers in the external momenta 
and/or meson masses. At this order many more terms contribute, especially counterterms 
which compensate their scale dependence. A complete calculation up to eighth chiral order, 
however, is far beyond the scope of this investigation and we will content ourselves by presenting 
numerical results at the scales fi = uip = 770 MeV and fi = m^/ = 958 MeV. The dependence 
of the non-analytic pieces on the scale might also give a hint on the importance of neglected 
LECs. Numerical results are given in Tab. |I|. 

Although the tree level contributions from the /9o,3,i3 term dominate the amplitude, the loop 
corrections turn out to be much smaller indicating the convergence of the chiral expansion. The 
same is true for the tadpole loops of fourth order couplings, they were found to be an order of 
magnitude smaller than the corresponding tree diagrams. A sample analytical expressions for 
the one-loop diagrams is relegated to App. y. 

4 Comparison with experiment 

Our numerical results for the decay widths must be compared with the experimental ones 

r(r/' ^ r/TT+TT") = (90 ± 8) keV 

T{7]' ^ r77r°7r°) = (42 ± 4) keV (34) 
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The pion and kaon masses are taken to be = 138 MeV, mj^ = 496 MeV and the pion decay 
constant is = 92.4 MeV. 

The large uncertainty in the LECs, in particular for /9o,3,i3) prevents us from making sensible 
numerical predictions. We have rather used the experimental decay width to pin down the LEG 
combination /5o,3,i3 more accurately. Taking the expression Aj^o up to second order ([l3|) one 
obtains /3o,3,i3 = (1-1±0.4) x 10"^ to be compared with the value /3o,3,i3 = (1-2±1.1) x 10"^ from 
conventional ChPT. The error bars could be reduced considerably while the central values are 
consistent. At fourth order in the chiral expansion of the decay amplitude loops enter together 
with new coupling constants. A fit to the decay width is possible but not unique. Further 
information is contained in the energy dependence of the decay which will provide a better test 
for our approach. To this end, two Dalitz variables x and y are introduced with 

X = ^{Et,i - E^2) 

2 + m^/m^ 2 + m^/m^ 
y = g (E^i + Kaj-lH g [m^i-mrj) (35) 

with Q = m^i — rrir, — 2mT, and -E^j being the energies of the pions. As Et,i and -E^2 vary 
over the physical region, x ranges from about —1.4 to 1.4 and y ranges from —1 to about 1.2. 
Experimentally, the decay of the rj' is then found to be described by the matrix element 

\A\^ = h{l + 2ay + ay"^ + cx^) (36) 



and fits to the data lead to a = -0.08 ± 0.03 for r]' t^tt+tt- |3| and a = -0.058 ± 0.013, 
a = ± 0.13^, c = 0.00 ± 0.03 for rj' rjii^Ti^ [^]. More recently the decay r]' —>■ rj-K^-K^ was 



used as a normalization for rare rj' decay searches and found to deliver a = —0.021 ±0.025 |33 
a value consistent with, but smaller than previous experiments. 

These experimental data must be compared with our calculations. We will keep the notation 
of (^). Since the /?o,3,i3 contribution dominates the amplitude, it is instructive to investigate 
its impact on the parameters in (p^). We do this by keeping only the /3o,3,i3 term in the 
amplitude while neglecting all other contributions. With these simplifications b in (^) is fixed 
by the value of /3o,3,i3 and the dependence of A on the Mandelstam variables determines the 
remaining parameters. The values of these slope parameters are functions of the masses of the 
four involved particles only and we obtain a = —0.138, a = —0.0016, c = —0.082 which turn 
out to be larger in magnitude than the experimental data. Albeit we are able to explain the 
measured decay width with the /3o,3,i3 term only, it does not reproduce the value for a. 

Let us therefore focus on the effect of further terms in Aj^o on the slope a. Most of those 
terms in App. ^ are constant in phase space and influence a indirectly since a is the relative 
slope. The only operators besides /3o,3,i3 yielding a sizeable slope are the and /^g^'' terms (the 
contribution from the 887r7r amplitude is suppressed due to mixing.) We find that the relative 



slope {a/ A in Tab. 0) of these terms is opposite and much larger than that of /5o3 13. Therefore 



these terms reduce a without affecting the decay width. By setting l3^^ = (3f^ = —0.7 x 10 ^, 



(3^^^ = /Sg^'' = 1.9 X 10 ^ and /?o,3,i3 = 1-1 x 10 ^ and all others at their central values our 
findings for the slope parameters resulting from are 

r = 42.3keV, a = -0.057, a = -0.022, c = -0.10. (37) 



■^The linear fit in Fig. 4 of [Q, however, shows a slope of a = —0.038, in contradiction with the value quoted 
and closer to the value of 1331. 
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Including the chiral corrections from Anlo ( |33| ) some of the LECs must be modified within 
their phenomenological ranges to restore the experimental decay width. Such a fit to the decay 
width is not unique. One possibility is to increase the value of /9o,3,i3 order to bring A^lo 
to agreement with the experimental decay width, while keeping the other couplings fixed. The 
slope parameters read then 

r = 43.6keV, a = -0.113, a = -0.019, c = -0.19 (38) 

at the scale /i = but one should keep in mind that a depends strongly on the choice 
of parameters and no precise statement about its value can be made at this level. E.g. by 
modifying some of the operators constant in phase space in order to accommodate the decay 
width instead of increasing /3o,3,i3 - the resulting \a\ tends to be smaller than the value shown 
in (HD. 

Under reasonable assumptions on the appearing LECs, such as neglecting OZI violating 
corrections to the couplings, we are thus able to achieve reasonable agreement with experi- 
mental data. Hence, (|37D constitutes our central values for the slope parameters and (^) is a 
conservative estimate of their uncertainties from A^^q related to the choice of the LECs. We 
observe a tendency to larger values of |a| if A^^q is taken into account. A non- vanishing slope 
I a I may represent the contribution of a gluon component to the Vj decay |3^. Fig. |^ shows the 
dependence of \A\^ on both, x and together with the experimental data from |Q and 

Furthermore, the leading order contributions to the slope a are due to counterterms of the 
fourth order chiral Lagrangian. In order to obtain a more reliable value for a, one must evaluate 
next-to-leading order chiral corrections including terms from the sixth order Lagrangian which 
we do not consider here. 

Nevertheless, it would be preferable to have an independent confirmation of the value for 
/3o- Such a check is provided by the calculation of scalar meson resonances using a coupled 



channel approach [^. In the present context, it is sufficient to sketch briefly some of the 
results, for details the reader is referred to [Q. In S'?7(3) ChPT the resonances /o(980) and 



ao(980) have been reproduced using coupled channels and taking only the second order La- 
grangian into account [^. In we applied a similar coupled channel analysis but using 



f7(3) ChPT and with the inclusion of the fourth order Lagrangian. While the established res- 
onances at 980 MeV remained unchanged, it gave rise to additional resonances in the isoscalar 
and isovector channels which could be identified with the /o(1370) or /o(1500) and ao(1450), 



respectively. This result is similar to where the fourth order Lagrangian had to be included 
to obtain the vector resonance p. Our results for the partial decay widths of the isovector agree 
with the experimental data of a^. The counterterms of the fourth order Lagrangian with four 
derivatives, such as the /Sq operator, turned out to be essential in producing the new resonances. 
In particular we found that the resonances were not seen for /5o,3,i3 < 0.5 x 10^'^ and did not 
fit the experiment for /3o,3,i3 > 1.5 x 10"'^. The other fourth order operators only changed the 
parameters slightly. This is actually not surprising since for momenta well beyond the scale of 
chiral symmetry breaking around 1 GeV momentum dependent terms will dominate (assuming 
they are comparable to other terms at smaller scales). In particular, the range for /?o,3,i3 is 
consistent with the value given in ( P^ and this confirms the choice we made for /3o,3,i3- 

Finally, we would like to comment on final state interactions due to two-pion rescattering 



as they have been evaluated, e.g., for the decay r\ — >• tttttt The final state interaction of 



two pions in the / = 0, s-wave channel is strong and attractive and is believed to dominate 
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Figure 2: Upper (lower) graph: The dependence of on the Dahtz variable y (x). The thick 
solid line denotes our result, and the dashed line represents the theoretical uncertainty due 
to the chiral corrections from A^i^o- Triangles and circles are the experimental data from 
and |33|, respectively. The thin lines are, descending in magnitude of slope, the linear fits from 
H and H (II). 
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interactions between the rj and a pion in rj' rjiT-K. In the Khuri-Treiman equations were 
used to determine numerically the unitarity corrections. The corrections to the decay amplitude 
were found to be 14% at the center of the decay region. We note that the energy of the two pion 
state is approximately the same for the decay rj' — > tjittt and roughly estimate the two-pion final 
state interactions for the present decay to be of similar size, i.e. less than 20%. Comparing this 
with, e.g., the uncertainty of the dominant coupling /3o,3,i3 (0) the corrections are moderate 
and smaller than the given uncertainty for /5o,3,i3. Hence final state interactions do not modify 
our conclusions considerably and we can safely neglect them. 



5 Summary and Conclusions 

In the present work, we have calculated the decay rj' rjTnr within the framework of infrared 
regularized U{3) ChPT. We presented the most general effective Lagrangian up to fourth or- 
der in the derivative expansion which describes the interactions at low energies between the 
members of the lowest lying pseudoscalar nonet, i.e. the octet of Goldstone bosons {TT,K,ri) 
and the corresponding singlet state r]'. At first sight, the fourth order Lagrangian looks less 
encouraging due to the proliferation of new unknown LECs. However, the unknown parameters 
entering the calculation can be identified either as OZI violating corrections of known LECs 
or couplings of parity-odd operators. Assuming that their contribution is negligible we are left 
with already known LECs and borrow their values from conventional SU (3) ChPT. 

The tree diagram contributions including 77-77' mixing have been evaluated first and it was 
found that together with the LEC combination v^^ from the second order Lagrangian, the two 
combinations /^s^ig and /5o,3,i3 from the fourth order Lagrangian constitute the decay amplitude 
up to second order in the Goldstone boson masses. In fact, it turns out that /?o,3,i3 dominates 
the decay amplitude and exceeds by an order of magnitude the contribution from The 
decay width can be used to reduce the uncertainty for /3o,3,i3 yielding /9o,3,i3 = (1.1 ±0.4) x 10~^ 
to be compared with /?o,3,i3 = (1-2 ± 1.1) x 10~^ from conventional ChPT. The dominance of 
/3o,3,i3 has been explained in the light of Adler zeros. Similar results at the tree level were 
obtained within large ChPT []T3| but loops were omitted since they are a next-to-next-to- 
leading order effect in the combined 1/A^c and chiral expansion. In contradistinction to large 
Nc ChPT we treat the 77' as a heavy field and evaluate loop diagrams by employing infrared 
regularization which keeps Lorentz and chiral invariance explicit at all stages of the calculation 
but suppresses contributions from loops with an 7/'. Within this scenario loop contributions 
start already at next-to-leading order. We restrict ourselves to the calculation of one-loop 
diagrams - both tadpoles and unitarity corrections - with vertices from C^"^^ and the LECs 
f^o'K f^A^^ and (3^\ In infrared regularization loop diagrams with a virtual rj' contribute only at 
higher orders and are therefore neglected. We find that the Goldstone boson loop contributions, 
in particular those with a /9o,3,i3 vertex, are rather small indicating the convergence of the chiral 
series. Reasonable agreement with experimental data is obtained and both the decay width 
and its energy dependence on the Dalitz variables are well reproduced. However, our results 
are subject to the choice of parameters such as the omission of OZI violating couplings. 

Overall we can state, that infrared regularized U{3) ChPT is suited to study the rj' decay 
7]' rjiiTi. We were able to reproduce the experimental data of the decay by taking into account 
next-to-leading order effects including one-loop corrections. Our results indicate furthermore 
the convergence of the chiral series. In the present work we did not discuss the possibility of 
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including low lying scalar mesons, especially the ao(980), as has been done in a tree level model, 
see [|16|. In our approach the explicit inclusion of these resonances is not necessary in order to 
obtain agreement with experiment although their effects might be hidden in some LECs. The 
ao(980) resonance, e.g., does indeed contribute to some of the couplings Pi in the resonance 
exchange picture [^. 
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A Abbreviations 

This is a summary of LEG combinations as they appear in the text: 



/5o,3,13 


oiO) I o{0) 3 niO) 
— Po ~r Ps ~ jPl3 


/34,5, 17,18 


= 3/?f ) + /?f - 9/?{? + 3P 


^5,18 






= - \V-Qv^' 








= /3f ) + IVgpS 


Pi" 








Pi" 





(0) 



B Amplitudes 

In this appendix, we present the amplitudes before applying rj-rj' mixing and r]' wave-function 
renormalization. As explained in the main text, the amplitude can then be divided into three 
pieces Aqstttt, ^887r7r and Aootttt with the subscripts denoting the external fields. The terms A2 
and A4 are the tree diagrams from the corresponding Lagrangians C^'^^ and C^^\ The terms 
Az, Amf stem from the wave- function renormalization of the octet fields and the expansion of 
meson masses and decay constants and A a, As^t,u are the tadpole and the unitarity corrections 
in the s, t, w-channel, respectively. The fourth order amplitude ^gsl-Tr been included in order 
to reduce the scale dependence, although it contributes via rj-rj' mixing to the rj' decay at sixth 
order in the derivative expansion. 
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3^ 

A4 + AZ + A^f + A^ + As + At + A^ 

J4 l"^"^7r "^r? j 3^ 

4\/2 



-{2ml, -s-t 



4x/2 



- I2{2ml + m^)/3f - 4(mJ + 2my /3f + + 

- 16(/3f ) - /3f ) (2m^ + m^) - 16(/3r - /jf 



45/.^ 



{2ml - 7ml) - 



2x/2m| 



,2~a) 



-(A^ + 5A,) 



8x/2m^4^)A, 2V2 4^)A^ 



+ 



+ 



— (2m, - 3s - m J 



2\/2 miv 2 I ■nnis) 2 2\ 

+ If, ' {^rnl-ml) 

2sf2 {2ml + mDvj'^AK AV2mlvi'h^,{t) 

9/6 

V2(2mUm>f)7..(t)^3^2^^2_3,^ 



+ 



2v/2 (2m^ + m^)i;^^^A^ 4v/2m;^i;^^^7^^(^) 



+ 



V2(2nil + n,^}r^-'l,,j^iu) 



(3m^ + m^ - 3m) . 
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Ai + Az + A^f 



A^ + As + At + A^ 



+ -?:r(^-"^r,-"^^')(«-2^J + IT 

Jtt '^Jtt Jtt 



4 MO) 



+ 



4^ 



3T 



J TT 



{t — ml — rn^,) {t — ml. — m^) + {u — m^ — m^,) {u — m^ — 

,-...2/3(0) 

m^(s — 2m^) + m^(s — — m^,) 



+ 



/4 



3(0) 



- 16(/3f - /3f - 16(4"^ - ) (2m^ + m^) - |A. - + 



?(o)^ 



?(o) 



(0)> 



+ 

ml 
ml 

^ (32m^ - 4m^ - 3s - 3t - 3«) - ^ (A, + 5 A,) 

9/4 36/4^ " ^ 
+ ^^^^ H - 2s) + ^^i^ {ml - Ami) 

+ ^^i^i^ + 3^ + 3mJ, - 9s) 
(3s + 18t + 3« - 12mJ) "''^^"^^^^ 



im /r- 2 I 2\ 



{2ml + ^l + ~ 3^) (3m^ + — 3t) 



36/4 
36/4 

^(3s + 3t + 18^-12m^)-^^%^ 
36/4 V 9/4 



Ikk{u) 
36/4 



(2m^ + ml + m^, — 3u) (3m^ + — 3ii) . 



,(2) Sv^^^ml 4-uf\ 2 2 \ 
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C A sample result for a higher order loop 

In this appendix, we present a sample result for the higher order loops which involve vertices 
from the Lagrangian C^'^^ We have chosen the contribution from a loop with two (3f^ vertices. 
The other contributions are obtained in a similar way but turn out to be rather lengthy. We 
have therefore omitted them for brevity. 



A.. 



) 6A^m^(52m^ + 221m2 + SOm^, - 33s) 



A, 



81/1 

+ 6AK{52m% + 19lm%ml - Ulm\ml - 18m^ 

+ 2>ml{2>m\ - 2m\){ml + 3m^) + 3s(9mt - llm^)) 
+ 2\ml{Sm\ - bml) (20m^ + ml + 3m5 + 3s) 

- l2h^{s)ml{Am\ + Sm^ + 3m2,) [l2,ml - 6s) 

+ Mkk{s) (2m^ + I2m\ml + 2ml - 3s(m^ + ml)) 

X (28m^ - I9m\ml + 3m2,(3m^ - 2ml) - ^s{m\ - ml)) 

+ AI^^{s)ml{%m\ - Sm^) [2m\ + ml) {Am\ - ml + mj,) 

^"^l^r^ [ - QKml{%ml + l^ml + 3m;, - 3t) 

- 6Ax(73m^m^ + 43m^mt + lOm^ 

+ 2>ml^,{m\ + 2ml){m\ + m^) + 9t{m]^ - Am\ml - m^)) 

- 6A,,mt(l6m^ + 5m^ + 3m;, + 3t) 

- 8V(t)mt(4m^ + 5m^ + 3m;,) (2m| + m^) 

- lKK{t){AQm]^ + 16m^m^ + lOm^ - 9t(5mi- + m^)) 



X 



{9m\ml + mj, (2m2 + m|) + 3t(m^ - ml)) (C.l) 
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